Laser-induced bound-state phases in high-order harmonic generation 
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We present single-molecule and macroscopic calculations showing that laser-induced Stark shifts 
contribute significantly to the phase of high-order harmonics from polar molecules. This is important 
for orbital tomography, where phases of field-free dipole matrix elements are needed in order to 
reconstruct molecular orbitals. We derive an analytical expression that allows the first-order Stark 
phase to be subtracted from experimental measurements. 
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It is a long-standing goal of atomic and molecular 
physics to follow electronic processes on their natural 
length and time scales, and eventually even to control 
their dynamics. Advances towards orbital tomography, 
the experimental reconstruction of an electronic wave 
function, have been made using high-order harmonic gen- 
eration (hhg) num. In HHG based tomography, a con- 
tinuum electron wave packet is formed at the peak of each 
half-cycle of an infrared driving laser. The continuum 
wave packet is accelerated in the oscillating laser field, 
and brought back to the vicinity of the bound-state wave 
packet at a later time. A high-energy photon is emitted 
when the continuum electron recombines into the ground 
state after having picked up kinetic energy in the laser 
field @, i]. The recombination step encodes informa- 
tion about the ground state onto the emitted harmonics 
through the recombination dipole matrix elements. 

In order to reconstruct the ground state orbital, one 
needs both the magnitudes and the phases of the recom- 
bination matrix elements. It is challenging to extract 
these phases from HHG spectra, as phases not related 
to the recombination step have to be subtracted from 
the harmonic phases. Phases related to the formation 
of a continuum wave packet and its subsequent acceler- 
ation in the laser field have already been accounted for, 
allowing the reconstruction of symmetric orbitals in N2 
and C0 2 0,(1. Nonsymmetric orbitals, however, pose 
a problem for orbital tomography due to their perma- 
nent dipoles, which cause them to acquire an additional 
laser-induced bound-state Stark phase. Using CO as a 
representative polar molecule, we find that the first-order 
Stark phase grows from zero to 0.57T within the harmonic 
plateau. The aim of this work is to extend the tomo- 
graphic method by accounting for this Stark phase, and 
presenting an analytical expression with which to sub- 
tract it from measurements. 

To isolate the effect of Stark phases we consider the 
idealized case of perfectly 3D oriented molecules. 

In the simplest model, the electric field F(i) of the driving 
pulse only influences the highest occupied molecular or- 
bital (HOMO) by adiabatically Stark-shifting its orbital 



energy E (atomic units are used throughout): 

E(F(t)) =E -»- F(i) - ±F T (t)^F(t). (1) 

As an example, we take CO oriented parallel to the po- 
larization axis of the driving laser. The HOMO has a 
symmetry, and an ionization potential of I p = \Eq\ = 
0.5150 au. Its permanent dipole = 1.1 au points 
from the carbon towards the oxygen nucleus. The compo- 
nents of the polarizability tensor a are a\\ = 3.2 au, and 
a±_ = 2.8 au [ll|. The resulting Stark shift caused by a 
moderately intense laser is sketched in Fig. [1] The chosen 
geometry suppresses contributions from the HOMO— 1 
due to its 7T symmetry [l^ . 

As the bound state evolves in time from t' to t, it accu- 
mulates the phase J* t * E(F(t"))dt". According to Eq. {!]) 
the accumulated phase differs from the field-free time- 
evolution by a Stark phase 
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where the first- and second-order Stark phases read 

^l Tk (t,t') = -j\-nt'')dt'' (3) 

*stlk(M') = -\J*F T (t")gF(t")dt". (4) 

The Stark phase changes the interference between the 
bound state, and the continuum wave packet that is born 
at time t' ', and returns at time t. 

Equations ©-((I]) are formulated in the time-domain, 
while orbital tomography is performed in the frequency 
domain. We therefore calculate classical electron trajec- 
tories, and translate the kinetic energy into photon en- 
ergy using energy conservation for those continuum elec- 
trons that recombine at the origin: 
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u(t,t') = -i 2 (t,t')+I p (F(t)). 



(5) 



Here r(t, t 1 ) is the velocity at time t of an electron that is 
detached at the origin with zero velocity at time t' , and 
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FIG. 1: (Color online) Sketch of the Stark shift felt by the 
HOMO of CO when subject to an electric field. The relative 
size of the first- and second-order Stark shift is drawn to scale 
for an electric field Fo = 0.075 au, corresponding to a laser 
intensity of 2 x 10 14 W/cm 2 . (a) The electric field is parallel 
to the permanent dipole of the orbital, and the effective ion- 
ization potential is raised, (b) The electric field is antiparallel 
to the permanent dipole, and the effective ionization potential 
is lowered. 



2p(F(t)) = \E(F(t))\. The effect of the ionic potential 
on the continuum trajectories is ignored. The result of a 
classical trajectory calculation is shown in Fig. [2] for an 
ultrashort pulse. The electric field is assumed to have a 
cosine squared envelope, and an 800 nm sine carrier wave. 
The full width at half maximum is 0.72 times one optical 
cycle T, corresponding to a full duration of 2T. The max- 
imum of the electric field is Fo = 0.071 au, corresponding 
to a peak intensity of 2.0 x 10 14 W/cm 2 of the envelope. 
Only one half-cycle contributes to the harmonic emission, 
which ensures that the continuum electron only recollides 
from one direction. Controlling the recollision direction 
is essential for tomography of nonsymmetric states [l3[ . 
Together, Eq. ([5]) and the calculation leading to Fig. [2] 
provide the required map between (t, t') and harmonic 
frequency. 

Returning to the Stark phases, Eq. (J3j> can be inter- 
preted in terms of the integral of the force felt by the 
continuum electron, showing that the first-order Stark 
phase is directly proportional to the return velocity 



(6) 



With orbital tomography in mind, one might be inter- 
ested in minimizing the Stark phase. However, Eq. © 
reveals that the first-order Stark phase cannot be min- 
imized by varying the driving pulse, as it only depends 




FIG. 2: (Color online) Classical trajectories for CO subject to 
an ultrashort laser pulse, indicated by the solid black curve. 
Instants of ionization are marked by circles, and instants of 
recombination by squares. The vertical offset of a trajectory, 
consisting of a circle and a square, indicates the photon en- 
ergy of the emitted harmonic, calculated using Eq. ([5]). Open 
(blue) markers indicate short trajectories, and full (red) mark- 
ers indicate long trajectories. 



on the return velocity, and the angle 9 between the in- 
ternuclear axis and the laser polarization axis. Insertion 
into Eq. ([5]) gives 

*aLk(*> =sg n (M ' *(*, *')) M |cos(0)| 



x y/2(u(t,t')-I p (F(t))), (7) 

where sgn(/x ■ r(t,t')) keeps track of the direction with 
which the returning electron probes the bound state. The 
time-dependence of the Stark-shifted ionization potential 
introduces a small difference between short and long tra- 
jectories. If the field-dependence of the ionization poten- 
tial is ignored, then Eq. simplifies to 



£LkM~±Mcos(0)y^ 



(8) 



The advantage of Eq. (JSJ is that it gives an analytical 
prediction of the first-order Stark phase directly in terms 
of the harmonic frequency, w, rather than in terms of 
electron trajectories through Lu(t, t'). We return to a dis- 
cussion of the accuracy of this result below. 

The second-order Stark phase cannot be expressed as 
simply in terms of the harmonic frequency. This is be- 
cause Eq. (|4|) turns out to depend explicitly on the driv- 
ing pulse, as well as giving qualitatively different result 
for short and long trajectories. The scaling of the second- 
order Stark phase with respect to laser parameters can 
be estimated by neglecting the pulse envelope, and inte- 
grating Eq. (j4]) from the peak of a half-cycle t' — t poa k 
up to t — tp Ca k + 2T/3. The cut-off harmonics are then 
found to accumulate a second-order Stark phase propor- 
tional to Fq/wq = 4wqC/ p , where luq is the frequency of 
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the driving laser, and U p the ponderomotive potential. 
According to the usual cutoff rule, w max = 3.17£/ p + I p , 
the importance of the second-order Stark phase can be 
reduced by increasing the wave length of the driving laser 
while keeping the harmonic cutoff fixed. 

We use the Lewenstein model [14J to test the validity of 
the trajectory calculations. The first-order Stark phase 
is extracted by subtracting harmonic phases calculated 
with and without inclusion of the first-order Stark shift 



in the Lewenstein model ll|,|l5|- Figure Eta) shows that 
short trajectories account for the first-order Stark phase 
within the limits set by the classical model, which does 
not account for harmonics below the ionization threshold, 
nor beyond the harmonic cutoff. The phase spikes coin- 
cide with minima in the harmonic plateau. The minima 
are caused by the interference between short and long 
trajectories, and each is associated with a sharp varia- 
tion in the harmonic phase. The exact position of each 
minimum changes when the Stark shift is included, which 
explains the phase spikes. 

The analytical model of Eq. (8) slightly overestimates 
the first-order Stark phase due to the use of the field- 
free ionization potential. The error is therefore ex- 
pected to increase with increasing intensity. Comparing 
Eq. ||HJ| to Lewenstein calculations at the 21st harmonic 
for three different intensities we observe an error of 8% 
(1.5 x 10 14 W/cm 2 ), 10% (2.0 x 10 14 W/cm 2 ), and 12% 
(2.5 x 10 14 W/cm 2 ). The interference features also move 
with intensity, but the square root behavior is the same 
at all three intensities. The close agreement between 
Eq. dHJ and the long trajectory calculation is due to the 
extremely narrow pulse envelope, which causes the long 
trajectories to recombine at low field strengths as shown 
in Fig. [2J 

The second-order Stark phase is plotted in Fig. [3jb) . 
The Lewenstein result is obtained by including the first 
and second-order Stark shifts in the Lewenstein model, 
calculating the harmonic phases, and then subtracting 
the phase obtained when only the first-order Stark shift 
is included. The calculation agrees qualitatively with 
the short trajectory prediction. The trajectory calcula- 
tion shows that the second-order Stark phase increases in 
magnitude when the electron spends longer time in the 
continuum. Selecting the short-trajectory contribution 
is thus an additional way of reducing the importance of 
second-order Stark phases. 

We would like to stress the point that the simple be- 
havior of the Stark phase is due to the fact that only 
one half-cycle, and mostly the short trajectories, con- 
tributes to the high-order harmonics. When several sets 
of trajectories contribute, the total phase is the result of 
a coherent sum of harmonic amplitudes, which can cause 
large modulations on top of the trend in Fig. [3] Exper- 
imentally, the dominating trajectory is selected through 
phase-matching by adjusting the position of the laser fo- 
cus relative to the nonlinear medium [la. 
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FIG. 3: (Color online) (a) Accumulated first-order Stark 
phase of CO subject to an ultrashort laser pulse (parame- 
ters are given in the text). The connected black dots indicate 
the phase obtained from a quantum calculation. The dashed 
(blue) curve is the phase predicted from short classical trajec- 
tories, and the dash-dotted (red) curve the phase from long 
trajectories. The smooth (green) curve is the analytical model 
of Eq. §8§ , with the phase defined to be zero for harmonics be- 
low 7 P . (b) Accumulated second-order Stark phase. Dashed 
and dash-dotted curves indicate short and long classical tra- 
jectories. 



In order to uncover the interplay between Stark phases 
and phase-matching, we now present results of the cou- 
pled solutions of the Maxwell wave equation (MWE) and 
the time-dependent Schrodinger equation (TDSE). We 
solve the MWE in the slowly evolving wave approxi- 
mation as described in [lf|, and at each plane in the 
propagation direction we solve the TDSE in the Lewen- 
stein model to calculate the time-dependent dipole mo- 
ment ll|. The nonlinear medium is a 5 mm jet of CO 
molecules. The CO molecules are assumed to be per- 
fectly oriented as in Fig. [2j The gas density is set to 



5 x 10 cm . We use the same driving pulse as above, 
except for adding a Gaussian focus with confocal parame- 
ter b = 2.0 cm. The focus is placed 0.70 cm in front of the 



middle of the jet. The peak intensity is 3.0 x 10 14 W/cm 2 , 
chosen so as to give a peak intensity of 2.0 x 10 14 W/cm 2 
in the middle of the medium. Our results are insensitive 
to ionization due to the very low target density. 

The Stark phase is calculated by propagating the 
MWE twice, with and without including the first- and 
second-order Stark effect in the single-molecule Lewen- 
stein calculations. In each case we transform to the far 
field, apply a spatial filter that selects predominantly the 
central, short trajectory contribution to the harmonics, 
and transform back to the near field. In an experiment 
this would correspond to having an aperture or a refocus- 
ing mirror in the far field. Then we subtract the phases 
obtained with and without including Stark shifts. An av- 
erage over the final spot on the detector screen is made by 
weighting the Stark phase $stark(w, r) at a given radius 
with the strength of the harmonic |Fhhg( w i r )\ 2 '- 



($StarkM) r 



Jq 00 $stark(w,r) |F H HG(^,' , )| 2 ^dr 
Jo 00 |F H HG(^,? , )| 2 '"dr 



(9) 



The result, shown in Fig. EJa), compares qualitatively 
with single-molecule predictions based on classical trajec- 
tories calculated for a peak intensity of 2.0 x 10 14 W/cm 2 . 
The phase oscillations stem from the intensity depen- 
dence of minima caused by interference between short 
and long trajectories. At a fixed intensity the minima are 
very sharp, giving rise to sharp variations in the extracted 
Stark phase as in Fig. [31 In macroscopic calculations the 
intensity falls off along, and perpendicular to, the prop- 
agation axis of the driving laser, causing the minima to 
smear out. The Lewenstein model has been known to 



overestimate the importance of the long trajectories [17], 
leading to an exaggeration of these interference oscilla- 
tions in our calculations. 

We also present results for perfectly aligned CO 
molecules. Spectra are calculated with and without the 
Stark effect by solving the MWE separately for opposite 
orientations, and adding the resulting harmonics coher- 
ently at the end of the gas [l8| . This procedure is valid in 
the limit where ionization-induced reshaping of the driv- 
ing pulse is neglible. The harmonics are then propagated 
to the far field, filtered, refocused, and the Stark phase 
extracted. Equation ((9J is used to calculate the radially 
averaged Stark phase shown in Fig. 0|b) . If the first-order 
Stark phase from opposite orientations had canceled out, 
then the total Stark phase would have been similar to 
that in Fig.[3jb). Instead, the total Stark phase is similar 
to that of fully oriented CO molecules in Fig. 0|a). The 
reason for this is that the Lewenstein model favors ion- 
ization when the electric field points from carbon to oxy- 
gen 11] . The ultrashort pulse only allows one half-cycle 



to contribute, thus increasing the relative contribution 
from the orientation shown in Fig. [2] The nonvanish- 
ing first-order Stark phase in Fig. U underlines that po- 
lar molecules behave differently from nonpolar molecules, 
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FIG. 4: (Color online) Stark phases from a gas jet of CO 
molecules (parameters are given in the text), (a) The CO 
molecules are oriented as in Fig. [2] The connected black dots 
are phases obtained from a full simulation including macro- 
scopic propagation. The dashed (blue) curve is the single- 
molecule phase predicted from short classical trajectories, and 
the dash-dotted (red) curve the phase from long trajectories, 
(b) The CO molecules are aligned parallel to the laser polar- 
ization. The black dots include both orientations of the target 
molecules. Dashed and dash-dotted curves indicate short and 
long classical trajectories. 



even if their head-to-tail symmetry is unbroken. 

In conclusion, we have investigated the role of laser- 
induced bound-state phases in HHG. We show that first- 
and second-order Stark shifts may both contribute sig- 
nificantly to the phase of harmonics generated from po- 
lar molecules. Such Stark phases must be accounted for 
if HHG is to be used for orbital tomography. We find 
a simple analytical expression for the first-order Stark 
phase, showing it to scale as the square root of harmonic 
frequency. No simple expression is found for the second- 
order Stark phase, but it can be minimized by increasing 
the laser wave length while keeping the ponderomotivc 
potential fixed. Propagation of the Maxwell wave equa- 
tion confirms that Stark phases survive phase-matching 
in the target gas. 
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